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Recent time-of-flight measurements on muon neutrinos in the OPERA neutrino oscillation exper- 
iment have found anomalously short times compared to the light travel-times, corresponding to a 
superluminal velocity, v — 1 = 2.37 ± 0.32 x 10~ 5 in units where c = 1. We show that cosmological 
bounds rule out an explanation involving a Lorentz invariant tachyonic neutrino. At the OPERA 
energy scale, nucleosynthesis constraints imply v — 1 < 0.86 x 10 -12 and the Cosmic Microwave 
Background observations imply v — 1 < 7.1 x 10~ 23 . The CMB limit on the velocity of a tachyon 
with an energy of 10 MeV is stronger than the SN1987A limit. Superluminal neutrinos that could 
be observed at particle accelerator energy scales would have to be associated with Lorentz symmetry 
violation. 
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I. INTRODUCTION 



Recent time-of-flight measurements on muon neutrinos in the OPERA neutrino oscillation experiment have found 
anomalously short flight-times compared to the light-travel times pQ. The anomaly corresponds to a supcrluminal 
velocity, v — 1 = 2.37 ± 0.32 x 10~ 5 in units where c = 1. The simplest type of particle with superluminal velocities 
would be a tachyon, a Lorentz invariant particle with imaginary mass, although this possibility conflicts with the 
observations of (anti)neutrinos associated with the supernova SN 1987A, which imply u-l<2x 10~ 9 for 10 MeV 
neutrinos [2HS]. We show that cosmological bounds support the supernova result and also rule out a tachyonic 
explanation of the neutrino timing anomaly. We find, at the OPERA energy scale, v — 1 < 0.86 x 10 -12 from 
nucleosynthesis constraints and v — 1 < 7.1 x 10~ 23 from the Cosmic Microwave Background (CMB) observations. 
The CMB limit on the velocity of a tachyon with an energy of 10 MeV is stronger than the SN1987A limit. 



II. TACHYONS IN AN EXPANDING UNIVERSE 



Our treatment of tachyonic particles in an expanding universe follows Refs. [7] and [8]. The tachyon is a particle 
with imaginary mass ifi and speed v > 1, whose energy and momentum in a local inertial frame are given by 

P = ^(v 2 -iy i/2 . (2) 

We also find it useful to express the speed as a function of the energy, 

^ 1/2 
E 2 



'(*0= + • ( 3 ) 



The tachyon speeds up if its energy is reduced. 

Consider such a particle in a spatially-flat Friedmann- Robertson- Walker universe with cosmological time t and scale 
factor a(t). If we use v l to denote the velocity is the co-moving frame, then the conserved linear momentum is 

Pi = fia- 2 v l . (4) 

The energy and momentum of the tachyon are related by 

Pl f-E 2 = ^ 2 . (5) 

The energy of a tachyon in an expanding universe is therefore given by 

The same result was obtained was obtained in Ref. [8j using the law of addition of velocities in an expanding universe. 
The important feature here is that the energy of a non-interacting tachyon falls faster than the usual 1 /a rate for 
masslcss particles. At some point the energy falls to zero and the tachyons disappear from the universe. This may be 
interpreted as a tachyon-antitachyon annihilation event [7]. 

The foregoing analysis only applies to free tachyons. A tachyon in thermal equilibrium with ordinary particles at 
temperature T would have energy E sa T, as usual. If the tachyons come out of equilibrium, then the constants in 
Eq. [6] are fixed by the energy at the relevant decoupling time. 



III. COSMOLOGICAL BOUNDS 



We shall start off the cosmological bounds with nucleosynthesis constraints (see e.g. [9] and [10] for a review), and 
focus on events between neutrino decoupling at a temperature T woa ^ and element formation at T nuc . An effective 
neutrino number N v is defined in terms of the energy densities of neutrinos p v and photons p 7 . For times preceding 
electron-positron annihilation, 



(7) 
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There is an extra factor of (4/11) 4 / 3 on the right-hand-side after electron-positron annihilation. In the standard 
model, with three lepton families, we have N v — 3. 

The cosmological expansion rate is sensitive to the number of neutrino species, and therefore both the neutrino 
decoupling time and the neutron/proton ratio at the time of helium formation are dependent on N v . By modelling 
the primordial element abundances and comparing them to observations it is possible to set limits on AN V = N v — 3. 
Typical examples of these limits are AN^ > -1.7 (95% CL) [10] and AN^ > -1.0 (95% CL) [9]. Note that these 
results assume AN^ is constant from decoupling to helium formation. 

Free tachyonic neutrinos, the energy is given by Eq. ^ with T oc 1/a, and we have E v = T wca k at neutrino 
decoupling to fix the constants, 



i2 rp2 \ I/ 2 



MT) = n U*— + 7/f- 1 ) ■ ( g ) 

V 1 weak r J 

Neutrino oscillation experiments imply that all three neutrino species would have a similar tachyonic mass |llj . The 
neutrino energy density p v oc T^E V and the photon energy density p 7 oc T 4 , so that the effective number of neutrino 
species defined by Eq. Q is 

N v (T) = 3(l + 7j !^-^) ' . (9) 

V J weak 1 J 

Note that N v = 3 at neutrino decoupling, then it decreases to zero and remains zero after the tachyon annihilation. 
Given limits on AN U at the nucleosynthesis time, then we can invert (|9| to set limits on the tachyon mass 

M 2 <^ T T" uc7 "° a 2 k |AiU (10) 

weak nuc 

Typical values would be T wca k = 0.8MeV for decoupling, T nuc = O.lMeV for light element formation and lAA^I < 2, 

fi < 0.12MeV. (11) 

The bounds on the tachyonic mass translate into bounds on the speed of the tachyon at a given energy scale via Eq. 
(§. If we where to compare with the OPERA results pQ, at an energy of 28GeV, 

v(28GeV) - 1 < 0.86 x 10~ 12 . (12) 

The speed reported by the OPERA announcement corresponds to v — 1 = 2.4 x 10~ 5 . 

CMB observations can be used to give information about the neutrino density of the universe at later times than 
the nucleosynthesis era. The acoustic oscillations in the CMB spectrum are sensitive to the free-streaming of the 
neutrinos prior to recombination time. The WMAP 7-year data analysis, assuming a ACDM model, places a limit 
N v > 2.7 (95% CL) [T^l US]- An earlier analysis of the WMAP 5-year data using independent input on the values of 
the expansion rate Hq and the parameter eg gave a similar limit N v > 2.55 (95% CL) [14] . 



Consider the latest limit with AN U > —0.3. We can replace T nuc in Eq. ( 10 ) by the temperature of the universe 
at matter-radiation equality when the density perturbations start to grow, T cq = 0.74eV. The limit on the tachyon 
mass becomes 

fi < 0.33eV. (13) 

The limits on the speed of the neutrinos at 28GeV become 

i;(28GeV) - 1 < 7.1 x 10~ 23 . (14) 

We also find that w(10MeV) — 1 < 5.4 x 10~ 16 , which is better than the supernova SN1987A limits on the neutrino 
velocity, although this only applies to the case of Lorentz invariant tachyons. With such small tachyonic masses, the 
three neutrino species could have a combination of tachyonic and real masses and still be consistent with neutrino 
oscillation experiments (which measure differences in /x 2 ), but this does not affect the bound. 



IV. CONCLUSION 



We conclude that the cosmological bounds rule out the possibility of a tachyonic neutrino over a wide range of 
tachyonic masses and the limits are better than those obtained from supernova 1987A. As a consequence, a GeV 
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energy-scale neutrino that travels appreciably faster than the speed of light would need to be associated with a 
violation of Lorentz symmetry. It would be possible to take specific models of Lorentz symmetry violation (see 
[T51 [To] for reviews) and repeat the cosmological analysis to obtain useful limits on these theories. 
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